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ABSTRACT 


Correlation techniques can be applied to 
infrared detection systems to improve the signal- 
to-noise ratio of the radiation detector. Two such 
correlation systems are compared and evaluated; 
(1) a single detector synchronous rectification 
system, and (2) a two-detector cross-correlation 
system. The two systems are interpreted in the 
context of search and tracking operations, and a 
discussion of the transfer of target information 
through these systems is presented. An experi- 
mental study verifying the theoretical analysis of 
the two-detector cross-correlation system has 


been carried out. 
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SYMBOLS 


Defined by Equation (4.25) 

True signal component contained in S 

True signal component contained in s 

Reference signal waveform 

Frequency; bandwidth limits 

Target signal frequency 

Banewidth of the noise, A v(t) 

Equivalent noise bandwidth (defined by Equation 2.20) 
rms value of AN(t) 

Noise component contained in S 

Power spectrum of Av(t) 

Maximum value of P(f) 

Amplitude-normalized power spectrum of Av(t) 
Improvement in signal-to-noise ratio (defined by Equation 4.11) 
Improvement in signal-to-noise ratio for unit bandwidth 
Unit-impulse respcnse function 

Correlator output signal 

Defined by Equation (4.8) 

Signal from radiation detector entering correlator 
Integration time 

Time 

Transfer function of integrator network 

A constant equal to 1.781 

Root-mean-square fractional error in Q(T) 

Dirac delta function 

Amplitude of sinusoidal target signal 

Amplitude of sinusoidal reference signal 

Defined on pp. 8, 11, and 12 with various subscripts 
rms value of A v(t) 

Noise component contained in s 


Time, as distinguished from t 
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I. INTRODUCTION 


In the past decade a number of technical reports have been 
written on the uses and applications of the correlation function. The 
foundation for this development was laid by Wienes! in 1930 with the 
introduction of the correlation function and its Fourier transform, the 
power spectrum. Applications have been made in such diverse fields 
as astronomy (Ref. 2, pp. 120-126), Biophysica and infrared physisee 

One of the most promising applications of thc correlation function 
is in the detection of signals in random noise. The initial papers of Lee 
and Wicener: (1950) and Lee, Cheatham, and Wieeneee (1950) present 
a discussion of the autocorrelation and cross-correlation methods for 
the detection of signals in noise. A somewhat amplified version of this 
work has been done by Fano’ (1951) and Daveapore: (1951). A mathe- 
matical analysis of the Dicke’? radiometer together with another system 
using two independent receivers has been carried out by Goldstein! ® 
(1955). Still more recently TacobsenT (1957) has exainined a multiple 
receiver system, laying primary emphasis upon obtaining directional 
information from the received signal. In addition to the above, the 
Massachusetts Institute of Technology (MIT) Research Laboratory of 
Electronics (RLE) has published several progress reports (e.g., Refs. 
12-14) dealing with correlators and their use in the detection of signals 
below the noise level. 

The relationship between the use of correlation techniques and 
optimum filter techniques is one of domain. The original work of 


North’ ” (1943) and Van Vleck and Middieron.< (1946) in the theory of 
eave cacriote refer to References 
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optimum filters is reviewed in a paper by Zadeh and Raggazzini 
(1952). Also appearing in this paper is a discussion of the relationship 
between correlation techniques and optimum filter techniques. Theo- 
retically, no claim to superiority can be made for either method, 
inasmuch as stationary random noise contains the same information 
whether described in the frequency domain (power spectrum) or time 
domain (autocorrelation function). However, from a practical point of 
view, the operation of a system in the time domain may be more feasible 
for many applications. Comments to this effect may be found in Ref. 6. 
The purpose of this report is two-fold. The first part is a theo- 
retical analysis of two detection systems, the synchronous rectification 
system and the two-detector cross-correlation system. The function of 
these systems is to improve the signal-to-noise ratio of radiation de- 
tectors. The analysis followed in the above references (viz. 5-11) is 
generalized; particular attention is paid to the effect which the form of 
the power spectrum has upon the improvement in the signal-to-noise 
ratio. Each system is then interpreted in the context of the current in- 
frared systems being used for search and tracking operations. The 
second part is an experimental confirmation of the two-detector cross- 
correlation system by means of a correlator developed at the MIT Naval 


Supersonic Laboratory (NSL). 
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Il. THE SYNCHRONOUS RECTIFICATION SYSTEM 


A block diagram of this system is shown in Figure 1. The 
signal and noise waveforms included in the diagram indicate the suc- 
cessive operations which improve the signal-to-noise ratio of the radi- 
ation detector. For convenience, the target signal is shown as a sinusoid. 
Radiation from a target -- after passing through the chopper -- falls on the 
detector and is sent to the preamplifier, appearing in the output as the wave- 
form indicated in Figure 1. A reference signal of the same frequency and 
phase as the target signal is produced by a steady light source and de- 
tector combination. The noise in the reference-signal waveform is 
assumed to be negligibly small. The two waveforms are subsequently 
multiplied together and then integrated; the result being that the target 
Signal will increase at a more rapid rate with time than the root-mean- 
square noise. The combined operations of multiplication and integration 
(or averaging) are called cross-correlation. 

The system indicated in Figure 1 may have several alternative 
forms, For instance, the multiplier can be replaced by a phase-sensitive 
detector. The function of the phase-sensitive detector is simply to reverse 
the sign of all waveforms during the negative phase of the reference wave- 
form, and pass all waveforms unchanged during the positive phase of the 
reference waveform. 

Before proceeding with the analysis a few definitions must be made. 
The output of the radiation detector before entering the correlator will be 


written as 


s(t) =e(t)+ Av(t) , (2.1) 


where e(t) and Av (t) are the true signal and noise components, respective- 


ly. The output of the integrator can then be written as 
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T 
S(t) = ( s(t) e' (t)dt ; (2.2) 


where e'(t) is the reference signal waveform entering the correlator. 


Expanding (2.2), it follows that 


T T 
S (t) = ie e(t)e'(t)dt + v Av(t) e'(t) dt (2.3) 


The first term of (2.3) is the signal component in the output of the 
integrator. 
ay 
E(T) = ( e(t) e'(t) dt (2.4) 
0 


The second term of (2.3) is the noise component in the output of the in- 
tegrator. 


An evaluation of the second term of (2.3) requires that some 
of the properties of the noise Av(t) be known. We shall assume that 


the cell noise, photon noise, and circuit noise comprising Av(t) have 


_ the following general properties: 








a) Stationarity, Av(t}?® = Av(t+tg for all ty 





2, - 
b) <Av(t) ap = <Av(tP>, 





c) Ergodicity, lim<Av(t}* > = Av(ttP 
T—-0o 


The brackets denote a time average over the period T, whereas the 
bar indicates an ensemble (statistical) average. No attempt to justi- 
fy the above properties will be made. The assumptions are not unduly 
restrictive, and the experimental data obtained on the basis of these 
assumptions substantially confirm their correctness. A discussion 
of the above properties may be found in Reference 18 together with 
appropriate source material. 

It is well known that the output of a linear network can be ex- 


pressed by the equation 
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00 
AN(t) = \ An(t)Ry (t — 7) dtr , (2.5) 


-O 


where An(r) is the input waveform to the linear network and 

Ry (t — rT) is the unit-impulse response function of the network. 
This equation can be verified by a simple application of the 
Fourier integral theorem. To find the statistical mean-square 


value of AN(t) we make the following operations: 


00 oo 
AN(t)? = , An(T)Ro(t — +) dr { An(r')Ro (t — 7!) dr! ; (2.6) 
~0 -00 


0 00 
AN(t)? = i‘ i An(7)An(7r’)Ry (t — T)Rg (t — 7') drdr! 3 (2.7) 
-0 % -0 


where T' is introduced to distinguish the two independent integra- 
tions. Taking the statistica) average of (2.7) and making use of 


property b) above, we have 


00 00 
AN(t)* = S ‘i An(7)An(7T') Ro (t — T)Ro (t — r') dr dr! . (2.8) 
-00 “ -00 


The stationarity property permits us to equate An(7T)An(r') with the 
autc covariance function U(r — T') of the noise An(r); that is, 


00 0 
AN(t)” = i ‘ W(T — T')Ro (t — T)Ryg (t — r')drdr' ; (2.9) 
-00 “ ~00 


The ergodic prope-ty implies the equality of the autocovariance and 
autocorrelation functions. Moreover, since the power spectrum of 
a random process is defined as the Fourier transform of the autocor- 


relation function, we may write 


ee) 
W(t — 7') -{ @(f) cos 2nf(r — r') df (2.10) 


-~ 00 


Equations (2.9) and (2.10) require only that the power spectrum of 


the noise and the unit-impulse response function be known in order 
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to obtain the mean-square noise in the output of the network. 
The autocovariance function of a stationary random process, 


e.g. u(t), is defined by the equation 


U(r) = u(t) u(t + 7) (2.11) 


where the bar denotes an ensemble average. Some of the properties of 


the autocovariance function include the following: 


W(T) = vu - 7) (2.12) 


and 


V(r) < w(O)= uw, (2.13) 


where u2 is the ensemble mean-square value of u(t). The autocor- 


relation function is defined by 


T 
¥(7) = lim ars 7 tthe +7) at (2.14) 


To _ 


In view of the ergodic property already assumed, the autocorrelation 
and autocovariance functions are equivalent. 


The unit-impulse response function of the integrator network of 
(2.3) is 


Tr 
Ro(t—n) = J 5 (t—7)dt , (2.15) 


where 6(t — 7) is the Dirac delta function. Equation (2.15) is very 


easily evaluated to give 
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' OS tet (2.16) 
Bo (=e »>T>TorT<dO 3 


Inserting (2.10) and (2.16) into (2.9) and carrying out the two inte- 


grals inv and T' gives 


0 
ar=s ( a Senet (2.17) 
“= 00 


The form of the power spectrum & (f) depends upon the power spectrums 
(or autocorrelation functions) of the signals entering the multiplier and 
the transfer function of the multiplier. 

Since the two functions in the second term of (2.3) are inde= 
pendent, the autocorrelation function of the product Av(t)e'(t) can be 
written as 

V(r) = P(r) Yt (7), (2.18) 
where wy (7) and w'(r) are the autocorrelation functions of Av(t) and e'({t), 
respectively. Using the fundamental theorem that the autocorrelation 
function and the power spectrum constitute a Fourier transform pair 
(see Eq.\(2.10)) the power spectrum associated with (rT) can be writtenas 


co co 
&(f) = t { P(f) cos 2nfr ai} y (r) cos 2nft dr 3 (2.19) 
= 00 = 


where P(f) is the power spectrum of Av(t). For convenience a noise 


equivalent bandwidth will be defined as follows: 


co 
Af*= , p(f) df ; (2.20) 
‘= 00 
where p(£) =— : (2.21) 
max 


In addition the mean-square value of Av(t) will be written as 


CO 
vos ‘ P(f) df = P Ai* 2. (2.22) 
) 0 max 


Inserting (2.19) into (2.17) and using the results of (2.20) and (2.22), 


the mean square noise in the output of the integrator becomes 
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= v2 
AN® =———=6{T) ,; (2223) 
2r 
where 
2 - (1— cos 2mfT),¢ 
‘ ‘ v'(7) cos 2nfr Spt cos 2nfrdf | dr ad eee (2.24) 
e(T) =2= = a 
( p(f) df 
-2 


Making the definition N= AN* , the root-mean-square noise becomes 


N == Ne it). (2.25) 


The general equation for the signal-to-noise ratio obtained 
by synchronous rectification can now be written. Taking the ratio 


of (2.4) to (2.25), we have 





T 
Em ‘ gS ties (2.26) 
N v le (tT) (T) 


This equation leads to the following conclusions. The signal-to- 
noise ratio obtained by the synchronous rectification system is 
dependent upon three functions: 

1) the form and frequency of the target signal 

2) the form of the reference signal 

3) the power spectrum of the noise entering 


the correlator 


A more detailed discussion of the effect of these three functions upon 
(2.26) will be given in Section IV. 

In applying the analysis to a system having real components, 
deviations from the ideal transfer functions of the multiplier and 
integrator must be included in (2.26), This means that the Fourier 
series coefficients of e(t) and e'(t), and the power spectrum p(f) 


will take on a somewhat modified form. 
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lI. THE TWO-DETECTOR CROSS-CORRELATION SYSTEM 


A block diagram of this system is shown in Figure 2. The sig- 
nal and noise waveforms included in the diagram again indicate the suc- 
cessive operations which improve the signal-to-noise ratio of the radi- 
ation detector. Radiation from a target passes through the chopper and 
falls on both of the detectors. The signals are then amplified and ap- 
pear in the outputs of the pre-amplifiers as the waveforms indicated in 
Figure 2. The signals, containing both trne signal and noise components, 
are then multiplied together and integrated. 

The analysis of this system follows the same pattern as the pre- 


ceding one. The two signals entering the multiplier will be written as 


S, (t) = e, (t+ Av, (t), (3.1) 
and 
8,(t) = e2(t) + Av, (t), (3.2) 


where e,(t) and Av (t) represent the true signal and noise components, 


respectively. The output of the integrator can then be written as 


re. 
S(t) = ‘ 8 (t)s2 (t) dt (3.3) 
0 
or 


T 
S(t) = ( [e, (t)e,(t) + e, (tlAv,(t) + e,(t)Av, (t) + Av, (t)Avz(t)}dt . (3.4) 
“o 


The first term of (3.4) is the signal component in the output of the 


integrator, 
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£ 
E(t) = { ce, (t)e, (t) dt . (3.5) 
0 


The last three terms of (3.4) are the noise components in the output 
of the integrator, 


The noise from these three terms can be written as 


AN(t) = AN, (t) + AN, (t) + AN; (t) ; (3.6) 


where AN, (t), AN, (t), and AN;(t) are defined by the equations 


eh 
AN, (t)= ( e,(t)Av, (t)dt (3.7) 
0 
T 
ANA( =f e(tdygirya (3.8) 
and : 
T 
AN; (t)= f Av, (t\Av,(t)dt —. (3.9) 
0 


Since the above three random noise variables are independent, it 
can easily be shown that 
AN? = AN? + AN? + ANY . (3.10) 


The first two terms in (3.10) have already been determined (see 
(2.3), (2.23), and (2.24)), The result is that 


2 


NP =54r 6,(T) (3.11) 
where 
N? = AN? ‘ v? = Av? ; 
and 
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oo oo = 
‘i Ww 2(r)cos 2nfr k p, (f)cos antrat] af cos 2nfT) 4¢ 
= 00 c 


O(T)= “T= So ee he) 
{ p, (f) df 
= 00 
Similarly, 
ie 
N3=s*7 62.(T) , (3.13) 
where 


00 
co { 2 1 (7)cos 2ufr [S pz (f)cos 2nfr a dr (i ecoremith df 
0,(T)= -o 00 -00 


u (3.14) 
( Pz (f) df 


=o 


The autocorrelation functions W;(7T) and W2(7T) are those corresponding 
to e,; (t) and e2(t); p, (f) and pz (f) are the amplitude normalized 
power spectrums corresponding to Av, (t) and Av, (t). 
The last term of (3.10) can be found by altering the form of (2.23) 
and (2.24). If the function w'(7) is redefined as the autocorrelation func- 


tion of 4v,(t), then w'(r) can be written as 


Y\(r)= ie P,(f) cos 2nfr df ; (3.15) 


- 0 


where P,(f) is the power spectrum of Av,(t). It therefore follows that 


se v2 p2 
Nyv= SF e(T) , (3.16) 


where 


00 eo oo . oo i 
{( [S p, (flcos 2nfr df ( p, (f)cos 2nfr a| cos2sfr sr} (= cos emt) af 
S65) 6S 


=o 


= 00 


0, (T) = -(3.17) 


oo ) 
( Pp, (f) df ( p, (f) df 
“eo “= 00 
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Using (3.11), (3.13), and (3.16) the root-mean-square noise in the out- 


put of the integration can be obtained: 


= 41% @s(T) t= 0,(T) + oz an! : (3.18) 
mY 2 VY 2 


Taking the ratio of (3.5) to (3.18), there follows that 





T 
( e, (t) e, (t) dt 
E _w2 *o ‘ (3.19) 
DE ECT) 
where 
1 1 
6(T)= 03 Em Og ties ceint) ° (3.20) 
This equation should be compared with (2.26). The signal-to-noise 


ratio obtained by cross-correlation is dependent upon the form and 
frequency of the target signals and the power spectrums of the noise 
entering the correlator. In the following section a detailed expres - 
sion for the ratio E/N will be discussed. As in the previous system, 
deviations from the ideal transfer functions of the multiplier and in- 
tegrator must be included in (3.19), if the analysis is to be inter- 


preted in terms of real system components. 
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IV. IMPROVEMENTS IN SIGNAL-TO-NOISE RATIO 


The resuits of Sections II and III provide a general description 
of the improvement in the signal-to-noise ratio which can be obtained 
with each system. A more detailed description requires that assump- 
tions be made concerning the form and frequency of the target and 
reference signals as well as the form of the power spectrum of the noise 
entering the correlator. In the analysis which follows the synchronous 
rectification system will be treated first. 

The true signal component of (2.1) we shall assume to be a sine 


wave of frequency fy and amplitude e, 

e(t) =e sin 2nfyt ‘ (4.1) 
The reference signal will be taken as 

e'(t) =e«' sin 2tfot ‘ (4.2) 


This equation permits the determination of the autocorrelation function 


Ww '(r) appearing in (2.24). From (2.14) it follows that 


wi'(r) = a cos2ufyt é (4.3) 


Furthermore, the integral appearing in (2.26) can now be evaluated. 





T ' : 
a _ €€ _ sin4 wf) T 
E = : e(t)e'(t) dt = 3 (1 =r ato T T . (4.4) 


TR 307 13 


CONFIDENTIAL 


CONFIDENTIAL 


For fy T >10(< 1% error), 


ee! 


5 T : (4.5) 


Ee 





The remaining function needed to evaluate (2.26) is the power 
spectrum of the noise entering the correlator. Some difficulty is 
encountered in selecting a power spectrum which is sufficiently 
general in form and yet reasonable to handle mathematically. 


Perhaps the most satisfactory one has the following form: 


eis lo forf,<f<f,, —f,< f< -f, (4.6) 
Oo , for—f, <f<f,, f>f,, {<—f, : 


where f, and f, determine the bandwidth Af of the noise. A diagram 
of (4.6) is given in Figure 3. 
On the basis of (4.3) and (4.6), the function 6(T) of (2.26) may 


be evaluated. The result is that 
eg pe 12 1 a 1 be 1 (4.7) 
2(f, — £,) =f ) 2 (fo + f,) *; f, ) Gr eS + f,) a, — f, ) 2 (f, — fy) és 
4 £98 2mT(fp —f,) , cos 2nT(f, —f cos 2nT(fp + f 
2(f) — f, ) 2 (fF, ae ) DATS +h) 
cos 2mT(f,p +f A : 
— ey + wT Si 2nT (fo + f, ) — wT Si 2nT (fo + f, ) 


+ wT Si 2nT(fyp —f,) + wT Si 2nT(f,— fy » } 


where the function Si(x) is defined as 


x _ 
six) = ( a du : (4.8) 
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| The above function converges very rapidly to w/2 as x — o; for x > 10 
Si(x) = 1/2 (within 6 %). An exhaustive treatment of the function can 
be found in Ref.19, pp. 1 - 9. 


The Equation (2.26) can now be written explicitly: 


Ret he! (1 — Sin Att =) (4.9) 
V2 ° NO(T) 

where 0(T) is given by (4.7). The signal-to-noise ratio obtained 

by synchronous rectification depends upon a complex function of the 

chopping frequency, bandwidth, and integration time. 

A close study of (4.9) reveals that a number of important trends 
can be determined, E/N approaches zero as T* for very small T, 
and varies as the square root of the integration time for very large T. 
Between the limiting values of T oscillations occur in a somewhat 
complex manner. 

The effect of the chopping frequency fy and bandwidth Af upon the 
improvement in the signal-to-noise ratio is also significant. The most 
satisfactory way to determine this effect is to discuss the form of the 
power spectrum of the multiplied noise (see (2.19)). The reason for 
this approach is that the integrator network following the multiplier 
operates as a low-pass filter. The transfer function of an integrator 
network integrating over a period of time T can be shown to have the 
following form: 


W,(f) = 527 {sin 2nfT + j(cos 2nfT — 1}. 
In terms of this equation (2.17) becomes 


co 
N =f af) |WOleae (4.10) 
co 


The function | w.(£)]? is shown in graphical form in Figure 4, 
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| It is apparent from the above two equations (and from (2.17) that 
the mean-square noise in the output of the integrator is much smaller 
for power spectrums having components predominantly at high frequency. 
We now proceed to a discussion of the power spectrum of the multiplied 
noise. 
When (2.19) is evaluated three cases are observed in the result: 
(a) fo — f, < f2 — fo < fo + £,; (b) fg — fp < fp — f£, < fp + fy; and (c) fy — f, < 
fo +t f£, < f, —fy. The power spectrums corresponding to these cases are 
shown in Figure 5 a,b, andc. These graphs show that low frequency 
components appear in the power spectrum of the multiplied noise that 
did not appear at all in the power spectrum of the cell noise (see (4.6)). 
Furthermore, no manipulation of the bandwidth or signal frequency fy 
can eliminate these low frequencies. The most that can be done is to 
minimize the contribution to the integral in (4.10) due to the low fre- 
quencies, and push as much of the power spectrum as possible into the 
high frequency range. To accomplish this fy should be set equal to f,; 
the chopping frequency fg should be as high as possible, since (from 
Figure 5 a,b) the high frequency components in the power spectrum are 
dependent upon its magnitude. The results of prescribing these condi- 
tions is that the value of the integral in (4.10) will be minimized. 
Having discussed qualitatively the effects of integration time and 
power spectrum upon the root mean-square correlated noise, (4.9) may 
now be represented graphically. However, let us first define a 


quantity Q(T) as follows: 


qr) = (EN, (4.11) 


(e /v) 


Q(T) represents the improvement in signal-to-noise ratio obtained for 


a given integration time and noise bandwidth. Written in terms of 
| Q(T), (4.9) becomes 


in 4nf, T Ele 
aT) =e (1-58 0 ) = (4.12) 
Jz Lato a) 
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Figure 6 shows a graph of (:.12) for Q(T) > 1 and for various values 


of the parameter Of. The signal frequency fg has been set equal to 


the high frequency limit f, of the band pass. The function Q(T), as 


shown in the graph, increases as the square root of the time-bandwidth 
product. It should be pointed out that the dependence of Q(T) upon the 
square root of the bandwidth Af is somewhat misleading, since the 
initial signal-to-noise ratio is inversely proportional to the square 
root of the bandwidth. 

A graph of an additional function Q*(T) is shown in Figure 7. 


This function represents the improvement in signal-to-noise ratio for 


noise of unit bandwidth (f, — f; = 1). The band has been moved over 


the frequency spectrum by assigning values to the center frequency 

£. (fy) + f,)/2. We may conclude from the figure that for improve- 
ments of magnitude Q(T) > 1, the system operates nearly independently 
of where the bandwidth Af is located in the frequency spectrum. In 
other words, the advantage obtained with high chopping frequencies is 
negligibly small. 

Proceeding to the two-detector cross-correlation system, we 
shall make the same assumptions with regard to the target signals and 
the power spectrums of the cell noise as in the previous system. The 
true signal components in (3.1) and (3.2) will be written as 






e2(t) =e, sin2mfygt 





The autocorrelation functions appearing in 6,(T) and 6,(T) (see 
(3.12) and (3.14) can be obtained from (2.14): 










2 
Wy(r) = =4 cos 2mfg T 







’ (4.15) 





£22 - 
5- cos 20fy T ; (4.16) 












W2 (tT) = 





Using the above two equations, (4.13) and (4.14), the integral in 
(3.19) can be evaluated to give 
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“o 


For fy T > 10 (< one per cent error), 


T 
sf @ (tle (t) dt = bez Ta (4.18) 


0 
The power spectrums of the cell noise, p, (f) and p,(f), are assumed 


to be identical, and of the same form as (4.6); that is, 


<{f< =f. 2S , 
Breen ie jfor sd, < F< = pat ees, avis) 


With the above assumptions the equation (3.19) can be carried out 


explicitly. The result is 








E t sin4nf,) T \ €, « T 
ie 1- jh el en , 4.20) 
N x ( Tio yr V2 rO(T) ( 
where 
@(T) = 65(T)+—L 0,(T) +L, (7), (4.21) 
v,? v2? 
and 


1 (f, + f,) 
63 (T) = —————_ & — 2lny — In4f,; f, + Zin ea Fy — 2ln2nT 
2(f2 — f, )? at 


+ cos 4nf, T + cos4nf, T + 4sin2nf, Tsin2nf,; T + Ci(4nf, T) 
+ Ci(4nf, T) — 2Ci2n(f, + f,)T +2Ci2n(f, — f,)T + 4uf, TSi4nf, T 


— 4nf, TSi2n(f, + £,)T + 4nf, TSi2n(f, — f, in} (4. 22) 
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eee 1 1 _ l _ i 
62(T) Tees ae, {ore Ei Cee = Cie ele tot 


4 COs 2uT(fyp — f) ) + COS 2nT(f, — fp ) 4 cos 2mT(fyo + £2) 
(4.23) 


eSOs Ca Ege ds or Sint iy eth wT Silat ley et) 
2(fo + f; ) 


+ nTSi2nT(fp —f,) + wT Si2nT(f, — a 


meet l = l - 1 a 1 
ONT) 2 ay att) 2 Pal Cig) el 


4 COS 2nT(f, — f£ es cos 2uT(f, — f 4 cos 2mT(fp + fz) 
iii (oh ee (4.24) 
— ns A) + wT Si2nT(fp + f,)— uTSi2nT(fy + f,) 

1 


The constant y is equal to 1.781. The function Ci(x) is defined as 
~~ cosu 
Ci(x) = — ( ——— du (4.25) 
~x u 


As x becomes large the function Ci(x) converges rapidly to zero as 
(sinx)/x; for x > 10, Ci(x) < 0.1. The function is treated fully in 
Reference 19, pp. 1 - 9. 

As in the previous system the rate at which E/N increases 
depends upon a complex function of the chopping frequency, band- 
width and integration time. However, some qualitative remarks may 


be made. For very small values of T the signal-to-noise ratio goes 
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to zero as T?; for large values of T, it increases as the square root 
of T. Between these limiting values no general characteristics can 
be obtained since the function behaves complexly. 

Let us consider the effect of the chopping frequency and the 
power spectrums of the cell noise upon (4.20). The three components 
of the noise in the output of the integrator stage are given in Equa- 
tions (3.6) through (3.9). The first two components in (3.6) are cross 
terms involving the target and cell noise signals. The power spectrums 
of each of these terms are nearly identical to those in Figure 5a,b,and 
c; the only changes are in the substitution ofe€, ore, for «' and P 


1max 
or P for eae These power spectrums are shown in Figures 
2 


8 and ‘9. The third component of (3.6) has a power spectrum as 
shown in Figure 10. The power spectrum of AN(t) is, of course, 
the sum of the three component spectrums. 

From the above three figures some conclusions can be reached 
for optimum correlator efficiency: 1) the chopping frequency should be 
kept very near the upper band-pass limit; that is, fy ~ f,;; and 2) the 
chopping frequency should be as high as possible. These conclusions 
follow from the fact that the integrator stage operates as a low-pass 
filter (see (2.17)). 

If the signal-to-noise ratios e, /v, and ez /v, are assumed to 
be equal (defined as ¢/v), then (4.20) may be written in terms of the 


function Q(T). The equation becomes 


aT) = (- Saat) & Saaz (4.26) 
where 0(T) is given by (4.21) through 4.24). The only significant dif- 
ference between this equation and (4.12) is that an additional factor 
e/v has appeared: Q(T) is now dependent upon the Original signal-to- 
noise ratio. If the ratio e/v is set equal to unity, then a graph of (4.26) 
may be made similar to that in Figure 6. The graph is shownin 
Figure 11. The signal frequency fy has again been given the same 
value as f2. As in Figure 6 the function Q(T) increases as the square 


root of the time-bandwidth product, for Q(T) > 1. However, it should 
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be remembered that the original signal-to-noise ratio e/v is in- 
versely proportional to the square root of the bandwidth. 

A graph of Q*(T) for the two-detector cross-correlation system 
is shown in Figure 12. The signal-to-noise ratio e/v has been set 
equal to unity; and fy, to f,;. As in synchronous rectification the improve- 
ment in signal-to-noise ratio for Q(T) > 1 is nearly independent of where 
the bandwidth is located in the frequency spectrum. That is, no signifi- 
cant advantage is gained by going to high chopping frequencies. 

Several conclusions can be obtained by comparing Figures 6 and 
ll. The curves in Figure 11 show a constant displacement downward 
of a factor of 2.45 from the curves in Figure 6. The cause of these 
lower Q(T) values are the additional two terms in the noise output of 
the correlator (see (3.6)).. A second conclusion relates to the dependence 
of Q(T) in two-detector cross-correlation upon the original signal-to- 
noise ratio. If another value for e/v were selected for Figure 11 th 
curves would be displaced. The trends in the displacement will be the 
following: 1) for e/v <<1,Q(T) is proportional to the first power of ¢/v; 
and 2) fore/v>>1,Q(T) is independent of the original signal-to-noise 
ratio. Thus synchronous rectification and two-detector cross-correla- 
tion are nearly equivalent fore/vy >>1. For e/v << 1, the improvement 
in signal-to-noise ratio in cross-correlation is less than that in syn- 


chronous rectification by a factor of 2/(e/v). 


TR 307 21 


CONFIDENTIAL ; 











CONFIDENTIAL 


V. RETICLE AND SCANNING SYSTEMS 


An important topic to be considered is the application of corre- 
lation techniques to reticle and scanning systems performing search 
and tracking operations. The primary function of these systems is to 
indicate the presence and position of the target within a specified field 
of view. This section discusses the transfer of positional information 
by the synchronous rectification and two-detecter cross-correlation 
systems. 

In reticle systems the position of the target is obtained from 
information inherent in the signal waveform. The sources of informa- 
tion are the phase, frequency, and amplitude of the signal. Various 
combinations of these sources are used to locate the target in the 
field of view. For instance, in the Sidewinder system the phase and 
amplitude of the signal provide the angular and radial information 
necessary for locating the target. When applying correlation techniques 
to reticle systems, the immediate problem is how to transfer positional 
information in the target signal and yet retain a relatively high informa- 
tion rate. Since positional information depends upon the shape of the 
signal waveform it necessarily follows that the correlation time must 
be short compared with the period of the signal. Otherwise, the origi- 
nal shape of the target signal is lost. The difficulty in obtaining signifi - 
cant improvements in signal-to-noise ratio for short correlation times 
( ~ 10°” sec) can be easily seen from Figures 6 and 11. Only for sys- 
tems having wide noise bandwidths (> 10kc) is this possible; for nar- 
rower bandwidths, the improvements are very small. Improvements 
in signal-to-noise ratio are accompanied by a reduction in information 


rate. We may, therefore, make the following conclusion: in reticle 
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systems requiring high rates of position information, it would appear 
that correlation techniques will find no application, unless the condition 
(Af/fp) 2 10° is met. 

The scanning systems referred to are those using the scan- 
ning spot and multiple-scanning spot. The single spot obtains target 
position information by systematically scanning and sampling small 

elements within the total field of view. The small element repre- 

| senting the instantaneous field of view localizes the target in the 
scanning pattern. The multiple-scanning spot is a linear array of 

single spots performing a one-dimensional scan to cover the field 
of view. Several detectors are used inthis system, each one cover- 
ing a portion of the scanning pattern. In both of these methods, the 
required improvement in signal-to-noise ratio and noise bandwidth 
of the system will determine the minimum dwell time of the scan- 
ning operation. The dwell time is to be equated with the correlation 
time. The design of the scanning system will then consist of a compro- 
mise between the parameters of information rate, total field of view, 
and instantaneous field of view. If low rates of target position informa- 
tion are acceptable in the design of a scanning system, then application 
of correlation techniques should give significant improvement in the 
signal-to-noise ratio of the system. 

The multiple scanning-spot method implies the use of more 
than one correlator in the scanning operation. This, in turn, raises 
the question of the size and weight of correlation systems. Early in | 
1960 work was begun on the producti of a miniaturized correlator. 
The result was a correlator composed of a Hall-effect multiplier and 
an RC filter acting as the integrator. A photograph of the system is 
shown in Figure 13. More than just the correlator is shown in the 
photograph. The multiplier unit is located in the upper right-hand 
corner; the integrator unit is in the lower right-hand corner. The 
elements on the left side of the photograph are transitor amplifiers 
which drive the Hall-effect multiplier. A free-running multivibrator 
circuit coupled with a relay was used to control the integration time 


over a period of about 0.25 seconds to 5 seconds. The scale at the 
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bottom of the picture is in inches; the third dimension of the correlator 

is about two inches. Figure 14 shows a circuit diagram of the correlator 
components and transistor amplifiers. The frequency response of the mul- 
tiplier is limited only by the characteristics of the coil, since the res- 
ponse of the Hall element itself extends well into the megacycle region. 

In summary we may conclude that correlators may be obtained which 
occupy a volume of four or five cubic inches and have a weight of a few 
ounces. 

One additional application of correlation techniques should be 
mentioned. A correlation system may be used to replace the cooling 
system of certain detectors; that is, detectors whose signal-producing 
mechanism is not impaired at room temperature. However, it should 
be emphasized that the use of a correlator will considerably increase 
the response time of the detection system. Recently, a successful ex- 
periment of this type was carried out with photoconductive InSb cells 
operating at room temperature. With a noise bandwidth of one kilocycle 
and an initial signal-to-noise ratio of unity, an order of magnitude im- 
provement in cell detectivity may be obtained in a correlation time of 
0.20 seconds. An order of magnitude improvement in detectivity is 
equivalent to cooling the cell to dry-ice temperature (195°K). See 
Reference 20, p. 1467 f. 


TR 307 25 


CONFIDENTIAL 


CONFIDENTIAL 


VI. THE EXPERIMENTAL STUDY 


Many correlators have been constructed for the purpose of uti- 
lizing the properties of the corzelation function. Among these may be 
distinguished two classes: the sampling correlator and the continuous 
correlator. The first class of correlators operates by sampling each 
time dependent function entering the multiplier stage over a given 
period of time; then the corresponding samples are multiplied and added 
together. Correlators of this type include those of Cowley”! (1949), 
Singleton”7(1950), and Cheatham”7(1951). The second class of corre- 
lators performs a continuous operation of multiplying the two functions 
together and then integrating (averaging) the result. Stevens~*(1950), 
Fano~(1950), Faran and Hills” (1952), and Johnson“ "(1952) have 
described correlators of this type. The correlator developed at NSL 
belongs to the second class. 

The two-detector cross-correlation system, incorporating the 
NSL correlator is shown in schematic form in Figure 15 and in pic- 
torial form in Figure 16. The radiation source is a Sylvania glow 
modulator tube, Type R1131C, modulated by a sinusoidal waveform. 
The modulated radiation source is a convenient component for the sys- 
tem since it eliminates the necessity of a chopper and provides a target 
signal of known form. The two PbS cells are mounted adjacently at the 
end of a brass tube. Since ac pickup is quite large in a high-gain network, 
the radiation source and detector assembly components are placed within 
a metal container (Figure 17). The preamplifier units following the PbS 
cells are Tektronix, Type 122, low-noise amplifiers. Bandwidth 
settings include the followirg: low frequency cutoffs (3 db), 0.2, 0.8, 8, 
and 80 cps; high frequency cutoffs (3 db), 50, 250, 1000, 10000, and 
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40000 cps. Gain settings are approximately 100 and 1000. The multi- 
plier stage of the correlator is a Philbrick, Model HM, combined mul- 
tiplier/operational fivefold. This hybrid component is designed for 
special purpose applications in analogue computing operations, and 
consists of a multiplier unit and five K2-W operational amplifier units. 
The multiplier has a dynamic range of — 50 to + 50 volts dc, a maxi- 
mum d-c error below one per cent of full scale, and a band pass of 
dc to 15 kilocycles. The integrator stage of the correlator uses one 
of the operational amplifier units, the remaining four being used as 
additional amplification stages beyond that of the preamplifiers. 

Oscilloscope pictures have been taken which illustrate ina 
graphical manner the operation of the correlator. Figure 18a shows 
a picture of the waveforms entering the two input channels of the cor- 
relator. The waveforms consist of a noise component alone: no radi- 
ation (target) signal is present. Figure 18b shows the output of the 
correlator, the integrator unit integrating successively over a period 
of 0.4 sec.; that is, the integrator is shorted out every 0.4 sec. The 
square divisions along the abscissa are calibrated to 0.2 sec per division. 
The zero line is located just below the central iine of the picture. Figure 
19a and b shows the change in the output of the correlator when a small 
sinusoidal target signal is introduced into the two input channels. Des- 
pite the small signal-to-noise ratio, the output of the correlator (Figure 
19b) shows in a remarkable way the presence of the signal component. 
The successive integrations are again carried out over a period of about 
0.4 sec. The scale along the abscissa is 0.2 sec per division, and the 
zero line is located three square divisions from the bottom of the pic- 
ture, 

To confirm the theoretical analysis of the two-detector cross- 
correlation system, (3.19) must be evaluated. Since the target signal 
is sinusoidal, the only functions needed to evaluate (3.19) are the ampli- 
tude-normalized power spectrums of the noise entering the correlator. 
These are obtained by means of a Briiel and Kjaer AF spectrometer 
(see Figure 20a,b). Since p, (f) and p,(f) are not analytic, a computer 


program was set up to evaluate (3.12), (3.14), and (3.17). The program 
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was carried out on an IBM 704 cornputer at the M.I.T. Computation 
Center. The results are shown in Figures 21 and 22 as the solid 
lines. The signal frequency fy was 500 cps, the initial signal-to- 
noise ratios were 1.0 and 0.50, and the equivalent bandwidths of the 
noise, Af,* and Af,”, were 1.16 kc and 1.04 kc. The signal-to-noise 
ratios entering both channels of the correlator were equal. 

From the experimental study, data points were obtained 
giving the improvement in the signal-to-noise ratio Q(T) in terms 
of the noise bandwidths, integration time, and initial signal-to-noise 
ratio. Two such series of data points are given in Figures 21 and 22. 

The low initial signal-to-noise ratios of 1.00 and 0.50 were 
produced by noting the linearity of the glow modulator, detector, and 
preamplifier system. The voltage across the glow modulator tube is 
proportional to the preamplifier output voltage over a wide range of 
voltages. In addition to the initial signal-to-noise ratio e/v, the two 
quantities required to determine Q(T) are E(T) and N(T). A statistical 
sampling procedure was followed in measuring these quantities, the 
samples being recorded on a Honeywell Visicorder. From the tapes, 
average and root-mean-square values of E(T) and N(T), respectively, 
were obtained. 

We proceed now to an error analysis of the experimental data 
points determining the function Q(T). Since Q(T) is defined by the 


equation 


Q(T) = EN) (6.1) 


it follows that 


AQT) JAE Av AN dc 
ci wae avenue ates tee) 


The root-mean-square fractional error in Q(T) is therefore 








See eee ee ee ae? 
- | 4Q7(7) AE? , Av? , AN? | Ae? 
An | —sSe2 = é 
eS |= * ae ae | i 
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The terms onthe right-hand side of (6.3) can be obtained by esti- 
mating the accuracy of the measurement of eachterm. The frac- 
tional error in measuring e and v was about 0.02 and 0.05, re- 
spectively. The number of samples taken in measuring E(T) and 
N(T) was 25. In the case of E(T) the error varies as a function of 
the correlation time, and has a range of values from about 0.05 to 
about 0.03. The fractional error in N(T) is approximately equal 

to the reciprocal cf the square root of the number of samples taken. 
With a sampling number of 25, the rms error is about 0.20. Thus 
the principle contribution to the error in Q(T) is due to the function 


N(T). Using the above estimates (6.3) can be evaluated to give 


A = 21 per cent (6.4) 


If the scatter of the data points around the zomputed value of 
Q(T) is assumed to follow a Gaussian probability distribution, 
then 68 per cent of the data points (or 2 out of 3) should lie within 
one standard deviation. Figures 21 and 22 agree with this assump- 
tion. 

We may conclude, therefore, that within the limits of ex- 
perimental error, substantial agreement was obtained between 


the theoretical analysis and the experimental study. 
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VII. SUMMARY 


In this report we have reviewed the possibility of applying 
correlation techniques to improve the signal-to-noise ratio of in- 
frared detection systems. General expressions have been derived 
which describe the operation of two correlation systems that can 
be used for this purpose. A detailed analysis for systems having 
band-limited white noise and sinusoidal target signals has also 
been presented. The most promising applications of these tech- 
niques appear to be with scanning systems and mosaics where 
high information rates are not required. An experimental study 
of the two-detector cross-correlation system has confirmed the 


theoretical analysis of this system. 
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Figure 3. 


Amplitude-Normalized Power Spectrum of Av(t) 
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Figure 4. Transfer Function of an Integrator Network 
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Figure 5. Power Spectrum after Multiplication (Synchronous 
Rectification) 
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Figure 8. Power Spectrum of AN, (t) 
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Figure 9. Power Spectrum of AN; (t) 
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Figure 10. Power Spectrum of AN; (t) 
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Figure 11. Improvement in Signal-to-Noise Ratio by Two-Detector Cross-Correlation 
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Figure 13. A Miniaturized Correlator Using a 
Hall-Effect Multiplier and RC Filter 


TR 307 47 


CONFIDENTIAL 





Siatjidury 10Ojstsuer], pue 10j;e[aIIOD sanzeuipy jo wrerseIg yMdIID “PT Vand 7 


4INDBID BOLVHD3ZINI 
? OS 











2 
aah 
<I 
}— wor 
z 
Li i701 HOZZ 
O 
LL 
=< t 
O ifeni 
i 
> Ag 
> Or 





elt ev 
3NNODUY 


48 


CONFIDENTIAL 








CONFIDENTIAL 


WIIOT DT}eUIaYdS UI UMOYS UTa\ysAg [eJUsUIIIadxXY OY ‘GT eansi 7 


cotter te rom KH 5 











Jarjpidwy 
HO1QIN, 


Jaryidwoasg 





XWOspMaL 








30199339Q Sag is 





hun 
Ko\dsig 


aayawousoeds ay 
sally Q janig 





| 


‘ 
10}2818G SQq fam 














wings sauyyduy JIISNOWOIs, 
Burjs0us HOUUQING x1U0sHOL 
30s0s6aj5uy [Ree ee eee OLE 


3030123305 


gounos 
volopoy 


PaOINpOW 





49 


TR 307 


CONFIDENTIAL 


WIO FT TEIIOJIq UI UMOYS uleyshg [eJUMUTTIEedxY ey ‘g{ erNn31qZ 


a 
= 
— 
rd 
LJ 
i 
Le 
Zz 
O 
O 


CONFIDENTIAL 





CONFIDENTIAL 


Ajquiassy 10499jaq pue xsoiNnoOg UOT eEIpey 


"LT eans1g 





51 


TR 307 


CONFIDENTIAL 


CONFIDENTIAL 





(a) Noise Components Entering Both Channels of the Correlator 





(b) Correlated Noise for Integration Time of about 0.4 sec. 


Figure 18 
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(b) Correlator Output Showing Presence of Signal. 


Figure 19 
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Figure 20. Power Spectrums of the Noise Entering Each 
Channel of the Correlator 
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Experimental Data Points and Theoretical Curve for 


(e/v) = 1.0 





Figure 21. 
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Experimental Data Points and Theoretical Curve for 


(e/v) = 0.5 


Figure 22. 
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